Surface metal-insulator transition in the Hubbard model 
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The correlation-driven metal-insulator (Mott) transition at a solid surface is studied within the Hub- 
bard model for a semi-infinite lattice by means of the dynamical mean-field theory. The transition 
takes place at a unique critical strength of the interaction. Depending on the surface geometry, the 
interaction strength and the wave vector, we find one-electron excitations in the coherent part of 
the surface-projected metallic spectrum which are confined to two dimensions. 
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I. INTRODUCTION 

The correlation-driven transition from a paramagnetic 
metal to a paramagnetic insulator is a fundamental 
problem in condensed-matter physics. The Mott transi- 
tion is of interest since strong electron correlations lead 
to low-energy electronic properties that are qualitatively 
different from those predicted by band theory. 

Theoretical attempts started with the early ideas of 
Mott 0, Hubbard @ and Brinkman and Rice f| (for an 
overview cf. Q). In recent years a more detailed under- 
standing has been achieved by studying the half-filled 
Hubbard model in the limit of high spatial dimensions 
d H . Here the electronic self-energy is a local quantity, 
and a self-consistent mapping onto the single-impurity 
Anderson model (SIAM) becomes possible [jsj Jj\. The 
impurity model can be solved numerically [0, |J |9|, |lO|, [ill] 
or, for weak or strong coupling, by using perturba- 
tive approaches (g, fL2[ . This constitutes the dynamical 
mean-field theory (DMFT) which is exact for d = oo but 
remains a powerful method also at finite d 13, fu) . 

Extensive DMFT studies for the d = oo Bethe lattice 
with a semi-elliptic free density of states have established 
the following scenario for the Mott transition |Q: For 
strong interaction U the two high-energy charge exci- 
tation peaks are well separated by an insulating gap in 
the one-electron spectrum similar to the Hubbard-III 
approach For decreasing U the insulating solution 
ceases to exist at a critical value U c \. On the other 
hand, for small U the system is a metallic Fermi liquid 
with a quasi-particle band at the Fermi energy. For in- 
creasing U the metallic solution of the DMFT equations 
coalesces with the insulating one at another critical in- 
teraction U C 2 > U c \. As U approaches U C 2 the effective 
mass diverges as in the Brinkman- Rice solution j3| . Be- 
tween U c \ and U C 2 there is a region where both solutions 
coexist. For T ^ the insulator is stable below U C 2, and 



the transition is of first order juj. For T = entropic 
effects favor the metallic solution. The transition is thus 
expected [|l4[ [ll| [ll| [lTj to take place at f/ c2 (second- 
order transition) or at least close to U C 2- 

Although some interesting questions are yet unsolved 
(e. g. concerning the concept of a pre-formed gap for 
U i > U C 2 Q), one can state that there is a compar- 
atively detailed understanding of the Mott transition 
in infinite dimensions. For finite dimensions one has 
to assume that the mean-field picture provided by the 
DMFT remains valid qualitatively. This is the basis for 
the present study which considers the Mott transition 
at T = 0, but concentrates on a new aspect of the prob- 
lem: the solid surface. Surface effects are of interest for 
different reasons: 

i) Taking into account the surface comes closer to the 
(inverse) photocmission experiment which determines 
the spectral function. The information depth in low- 
energy electron spectroscopy is usually limited to a few 
surface layers only JU^. 

ii) The breakdown of translational symmetry due to 
the surface introduces a layer dependence in the physical 
quantities which is worth studying. The layer depen- 
dence of the quasi-particle weight is calculated within 
the DMFT and also within conventional second-order 
perturbation theory around the Hartree-Fock solution 
(SOPT-HF) (20). The perturbative treatment applies to 
the weak-coupling regime, non-local terms can be stud- 
ied here. The (non-perturbative) DMFT accounts for 
the local (temporal) fluctuations exactly, but neglects 
all spatial correlations. Within both approaches it is 
found that the presence of the surface gives rise to strong 
oscillations depending on U and geometry. 

iii) One may ask whether or not the critical interac- 
tion strength is modified at the surface, analogously to a 
(possibly) enhanced Curie temperature at the surface of 
a ferromagnet plf. For uniform model parameters we 
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find a unique critical interaction U C 2 where the whole 
system undergoes the transition. This excludes a sur- 
face phase different from the bulk. 

iv) The low-energy electronic structure is of par- 
ticular interest close to the critical point. Here but 
also well below U C 2, the surface renormalization of the 
one-electron excitation energies is shown to be suffi- 
ciently strong to generate a surface mode. The mode 
is essentially confined to two dimensions and shows up 
with a reduced dispersion. The generating mechanism 
for this type of surface excitation is exclusively due 
to electron correlations. A simple criterion for their 
existence is derived analytically. 



II. MODEL AND GREEN FUNCTION 



We study the Hubbard model at half-filling and T = 0: 



(1) 



i j a 



with i and j running over the sites of a semi-infinite 
d = 3 lattice which is chosen here to be the simple cubic 
lattice cut at a low-index lattice plane. The uniform 
nearest-neighbor hopping tuj\ = —t with t = 1 fixes 
the energy scale. To focus on the Mott transition from 
a paramagnetic metal to a paramagnetic insulator, we 
restrict ourselves to the spin-symmetric solutions of the 
DMFT equations as usual 0] . Thereby, we set aside the 
fact that on a bipartite lattice the half-filled Hubbard 
model is actually in an antiferromagnetic phase below a 
finite Neel temperature. 

We concentrate on the metallic spectrum in the vicin- 
ity of the Fermi energy. This is governed by the low- 
energy expansion of the self-energy fl2^] : 
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With the help of the (symmetric) matrices (3 = (/%) and 
Z = (1 — /3) _1 the one-electron Green function Gij(E) = 

((cia-', c j a )) can b e written in the form: 



G(E) = Z T 



El - ZiTZi 
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where G = (Gy) and T = (ty). Only the term linear 
in E is taken into account. In particular, we consider 
an energy range where damping effects (ex ijijE 2 ) are 
unimportant. 

This coherent part of the spectrum is exclusively 
determined by the so-called quasi-particle weight matrix 



Z = (zij) (mass-enhancement matrix Z _1 ). The usual 
bulk quasi-particle weight z(q) is obtained from z^j 
by three-dimensional Fourier transformation. For the 
semi-infinite system only lateral translational symmetry 
can be exploited. Considering the system to be built up 
by N± layers parallel to the surface (N± i— > oo), there 
is a collection of N± one-dimensional Fermi- "surfaces" 
in the two-dimensional surface Brillouin zone (SBZ). 
The /x-th Fermi surface is given by — where ek M 
(fi = 1, ...,N±) are the eigenvalues of the renormalized 
hopping matrix ZaTZi at a two-dimensional wave 
vector k of the SBZ. The corresponding eigenvectors 
Uk M yield the discontinuous drops of the momentum 
distribution function along the direction k in the SBZ 
via: Sn(k = k£) = Nj x uL • Z • u kA1 . For N±_ i-> oo 
the decrease of n(k) is continuous, the quasi-particle 
weight matrix Z, however, has a well-defined meaning 
according to Eq. (ph. 



III. PERTURBATIONAL APPROACH 

Z shall be calculated at a mean-field level assuming 
that spin and charge fluctuations are reasonably local. 
This implies a local self-energy, Ey « Sij , which needs 
justification for d < oo and in particular for the d = 3 
semi-infinite lattice. Some insight can be gained by con- 
ventional second-order perturbation theory around the 
Hartree-Fock solution. SOPT-HF is capable of account- 
ing for the complete non-locality of the self-energy in 
the weak-coupling regime also for the case of reduced 
translational symmetry p(i| . 

We have calculated the linear coefficient for the 
different low-index surfaces. The low-energy expansion 
of the SOPT-HF self-energy yields: 
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(4) 

if the HF (on- or off-site) density of states is symmetric, 
Pij\~E) — Pij\E), and /3y = if it is antisymmet- 
ric. Fig. 1 shows the local term to exhibit a weak 
layer dependence for the sc(100) surface. Only for the 
topmost surface layer a significant difference to the bulk 
value is found. Surface effects become stronger for the 
sc(110) and are most pronounced for the sc(lll) sur- 
face where a strong oscillation of (3u is noticed. This 
can be understood as follows: According to Eq. (Q) the 
dominant contribution to (3n comes from the HF density 
of states density near E = 0. For the sc(lll) surface, 
however, the E = value is known p3| to oscillate be- 
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tween zero (for the sub-surface and all even layers) and 
a constant value (for the topmost and all odd layers). 
This anomalous behavior is found (23) to dominate the 
shape of the density of states within an energy range 
AE ~ tL~ 2 , where L denotes the distance to the sur- 
face. Passing from the surface to the volume, L *—>■ oo, 
the bulk HF density of states is recovered outside an in- 
finitesimally small energy range around E = 0. For the 
coefficient j3u this implies a damped layer-by-layer oscil- 
lation with the maximum absolute value for the topmost 
layer. 

For the discussion of the non-local terms ft^j we have 
to consider second nearest-neighbors since electron-hole 
symmetry requires p\j (E) to be antisymmetric and thus 
flij = for nearest neighbors i,j on a bipartite lattice 
and at half-filling. We also find surface-induced oscilla- 
tions in the case i ^ j as can be seen in Fig. 1. More 
important, however, the absolute values are smaller by 
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Fig. 1. Layer dependence of the on-site (left) and the off- 
site (right) linear coefficient in the low-energy expansion (^) 
of the self-energy at different surfaces of the s.c. lattice as 
obtained within SOPT-HF. Second-nearest neighbors i and j 
within the same layer. Within SOPT-HF /3y ~ U 2 x const.; 
we set U — 1 in the figure. Note the different scales. Energy 
units such that t = 1. 




Fig. 2. U dependence of the layer-dependent quasi-particle 
weight z a = (1 - Pa)' 1 within SOPT-HF. Results for differ- 
ent surfaces. '1' stands for the topmost surface layer. Dashed 
lines: bulk s.c. lattice. 

about two orders of magnitude compared with the local 
terms. 

Let us consider the quasi-particle weight matrix. De- 
composing the linear coefficient into local and non- 



local parts, Pij 
Z=(l-/5) 



Sijflij + (1 — and expanding 



1 in powers of the non-local one, yields 



Az u = 



1 



ft. 



1 



(5) 



for the lowest-order non-local correction Azu of the local 
element za. The correction is of the order (/3j^-) 2 and 
can thus be neglected in the weak-coupling limit. 

Fig. 2 shows the local elements of the quasi-particle 
weight matrix za sw 1/(1 — f3a) within SOPT-HF as 
a function of U. For weak coupling we observe the 
quadratic dependence on the interaction strength: 
1 — z a (U) oc U 2 . Depending on geometry, there is a 
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considerable layer dependence even for weak interaction. 
For strong U, SOPT-HF still predicts a Fermi-liquid 
state with z a > 0. Here, as expected, the perturbational 
approach breaks down. 

IV. DYNAMICAL MEAN-FIELD THEORY 

We now turn to the dynamical mean-field theory to 
consider the intermediate- and strong-coupling regime. 
DMFT is a non-perturbative method. However, we have 
to assume that a local self-energy functional is a reason- 
able approximation for d = 3 dimensions. From the 
perturbational results, this is well justified in the weak- 
coupling regime, for larger U the assumption may be 
questioned. In any case, we expect DMFT to be a good 
starting point. 

We have generalized the mean-field equations to the 
case of film geometries. A film consisting of a sufficiently 
large but finite number of layers N± is considered to 
simulate the actual surface. Apart from mirror symme- 
try at the film center, the layers are treated as being 
inequivalent. DMFT must thus be applied in the fol- 
lowing way: We start with a guess for the (local) self- 
energy H a (E) for each layer (translational symmetry is 
assumed within the layers). Second, the Dyson equation 
is solved on the semi-infinite lattice to get the on-site 
Green functions G a (E). For this purpose we make use 
of lateral translational symmetry which restricts the nu- 
merical calculations to inversion of N± x N± matrices 
on a wave- vector mesh in the two-dimensional SBZ. Cor- 
responding to each layer a, an effective impurity prob- 
lem is set up. From the on-site Green function and the 
self-energy, the a-th SIAM is specified by fixing the re- 
spective free impurity Green function G^ p (E) via the 
DMFT self-consistency condition: 

GiZ(E) = (G a (E)- 1 + Z a (E))- 1 . (6) 

The crucial step is the solution of the (a-th) SIAM to 
get the self-energy E Q (E) which is required for the next 
cycle. Applying DMFT to a film geometry implies a self- 
consistent treatment of N± impurity problems that are 
coupled indirectly via the respective baths of conduction 
electrons. 

For finite temperatures the impurity problems can be 
solved by employing the Quantum-Monte-Carlo method 
0, HH using the Hirsch-Fye algorithm. For T = the 
Exact-Diagonalization (ED) approach of Caffarel and 
Krauth [jio], |l) may be applied and will be chosen for 
the present study. ED is able to yield the essentially ex- 
act solution of the mean-field equations in a parameter 



range where the errors introduced by the finite system 
size are unimportant. For the Mott problem the relevant 
low-energy scale is set by the width of the quasi-particle 
peak in the metallic solution. It has to be expected that 
there are non-negligible finite-size effects when this en- 
ergy scale becomes comparatively small. We are thus 
limited to interactions strengths that are not too close 
to U C 2 and cannot access the very critical regime. 

The algorithm proceeds as follows: The DMFT self- 
consistency condition yields a bath Green function 
G^ p (iE n ). The parameters of a SIAM with n s sites, 
the conduction-band energies Ck and the hybridization 
strengths Vk (k = 2, ...,n s ), are obtained by minimizing 
the following cost function: 

1 "max 

'''max i 1 I 
n— 

(7) 

where G { £ (iE n - fi) = (iE n -e d -J2 k V*/(iE n - e,,))- 1 
is the free (U = 0) Green function of the n s -site SIAM. 
The fit of (iEn) -1 is performed on the imaginary 
axis, iE n — i(2n + l)ir/[3 with a fictitious inverse 
temperature (3 which introduces a low-energy cutoff. 
Lanczo's technique J24| is used to calculate the zero- 
temperature impurity Green function Gi mp (i-E n ). The 
local self-energy of the a-th layer is then obtained via 
Z a (iE n ) = G^iEn)- 1 - GWpfiE,,)- 1 . At half-filling 
electron- hole symmetry can be enforced by e& = — e^' 
and V£ = V£ with k + k' = n s + 2. 

V. RESULTS AND DISCUSSION 

Routinely, the calculations have been performed for 
n s = 8 sites. For interaction strengths well below U C 2 
this has turned out to be sufficient for convergence. 
The results are independent of the high-energy cut- 
off n max . A small low-energy cutoff (/3W ~ 10 3 ) is 
necessary to obtain a converged value for 1 — l/z a = 
(d / d(iE))lmT, a (iE = 0). Although the film geometry 
implies a high-dimensional parameter space (N±(n s — 
l)/2, i. e. ps 100 parameters to be determined self- 
consistently) we always obtained a stabilized metallic or 
insulating solution. A moderate number of layers N±_ in 
the film is sufficient to simulate the semi-infinite system. 
We used N± = 11, 15,21 layers to investigate the (100), 
(110) and (HI) surfaces, respectively. The convergence 
has been checked by comparing the results from calcu- 
lations for different N±. 

Difficulties arise for U i— > U C 2- Due to the finite 
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Fig. 3. Bulk quasi-particle weight as a function of U. Results 
for n s — 8 and n s = 10. 



number of sites considered, the energy resolution of the 
conduction band cannot be better than AE w rjW/n s , 
where W is the free band width and 77 is a constant which 
accounts for the fact that the conduction-band energies 
6k are not equally spaced and depend on U. AE can 
be estimated by comparing the results for different n s . 
Fig. 3 shows the bulk quasi-particle weight z as a func- 
tion of U as obtained from calculations for n s = 8 and 
n s = 10 sites. We notice that there is a good agree- 
ment between both results for U smaller than ~ 14.5 
or z larger than ~ 0.01. For z < 0.01 the energy scale 
set by the width of the quasi-particle peak w zW can 
no longer be resolved, i. e. AE w 0.12 = W/100. For 
n s = 8 and U = 15 this AE is also the typical ener- 
getic distance between the conduction-band level at the 
Fermi energy, et — 0, and the level with lowest absolute 
energy ey ^ in the self-consistent solution. JThe error 
that is introduced by the low-energy cutoff f3 is found 
to be of minor importance compared with the error due 
to finite-size effects. This is ensured by choosing (3 such 
that tt/?- 1 (w 0.03) is smaller than AE. 

From the zero of z(U) in Fig. 3 we can estimate the 
critical interaction: U C 2 w lQt — 1.33W for n s = 8 and 
U C 2 « 15t = 1.25W for n s — 10. Obviously, a pre- 
cise determination is not possible. The results may be 
compared with the value U C 2 ~ 1.34W which has been 
obtained for the Bethe lattice with infinite connectivity 
using ED at a finite but low temperature (3t — 100 p5[ . 
For interaction strengths U > U c \ = 11.54 we also ob- 
tain an insulating solution. Contrary to U C 2, the value 
found for U c \ is almost independent of n s . This is plau- 



sible since there is no small energy scale in the case of 
the insulator. 

We now return to our main subject of interest: the 
layer and interaction dependence of the quasi-particle 
weight at the surfaces of the s.c. lattice. Fig. 4 shows the 
results obtained for n s = 8. At weak coupling we notice 
a quadratic U dependence consistent with perturbation 
theory. For the sc(100) surface there are no significant 
differences between the DMFT and the SOPT result for 
z a up to U «2-3 (see top panel in Fig. 4). The same 
holds for the sc(110) surface (not shown). Contrary, 
perturbation theory breaks down at a much weaker in- 
teraction strength for the sc(lll) surface as it is obvious 
comparing the results in Figs. 2 and 3. 

For low and intermediate U, z a has an oscillating layer 
dependence while it is monotonous close to the transi- 
tion (Fig. 4, insets). The strongest surface effects are 
seen for the open sc(lll) surface while for the sc(100) 
surface z a ps ztmik except for a = 1. This trend is re- 
lated to the decrease of the surface coordination num- 
bers: rig 100 ' = 5, n.g 110 '' = 4, ".g 111 ^ = 3, to be compared 
with the bulk value rtB = 6. 

For all surfaces and layers there is a unique and com- 
mon critical interaction U C 2 where z a (U) approaches 
zero. The value of U C 2 is the same as the bulk critical 
interaction strength (found for n s — 8), i. e. U C 2 is deter- 
mined by the bulk system. In all cases the quasi-particle 
weight of the topmost layer zs is significantly reduced 
with respect to the bulk value z. This is plausible since 
the variance A = n$t 2 of the free (Z7 = 0) surface density 
of states is smaller due to the reduced surface coordina- 
tion number ns , and thus U / VA is larger at the surface 
which enhances correlation effects. The reduced surface 
coordination number therefore tends to drive the surface 
to an insulating phase at an interaction strength lower 
than the bulk critical interaction U C 2 (z > 0, zs = 0). A 
real surface transition, however, is not found; zs remains 
to be non-zero up to U C 2- 

This is interpreted as follows: For U < U C 2 the bulk 
quasi-particle band has a finite dispersion. The corre- 
sponding low-energy excitations are thus extended over 
the entire lattice. Due to hopping processes between the 
very surface and the bulk, they will have a finite weight 
also in the top layer. Therefore, below U C 2 the bulk exci- 
tations will to some extent induce a quasi-particle peak 
with a non-zero weight z$ > in the topmost layer. To 
test this interpretation, we decoupled the top layer from 
the rest system by switching off the hopping between the 
top and the subsurface layer, t 12 — 0. For the sc(lll) 
surface this implies n§ = 0. The top-layer self-energy 
is thus given by the (insulating) atomic-limit expression 
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Fig. 4. U and layer dependence of z a = (1 — /3a) _ in the 
vicinity of the surface (bulk: dashed lines). ED calculation 
for 8 sites (filled/open circles: n s — 10). N± = 11, 15,21 for 
the sc(100), sc(110), sc(lll) surface, respectively. Insets: z a 
in the critical regime. Thin (solid and dashed) lines: SOPT- 
HF results for the sc(100) surface and U < 8 (top panel). 



= U/2 + U 2 /4E, while for the other layers we 
get the same layer- dependent (metallic) self-energy as 
before but shifted by one layer. Taking this as the start- 
ing point for the DMFT self-consistency cycle where 
ti2 = t again, we observe that a finite quasi-particle 
weight in the top layer zg > is generated immediately 
and that the cycle converges to the same solution as 
before. 

There is no Mott localization at the surface of the 
metallic Fermi liquid. However, the tendency of the 
system to reduce the surface quasi-particle weight (en- 



hance the surface effective mass) results in a partial 
(wave-vector-dependent) localization of quasi-particles 
at the surface: Let us consider the coherent part of the 
surface electronic structure. To make the calculations 
more transparent, we assume that z a = z except for 
the top layer {z a =\ — z$). This is well justified for the 
sc(100) and, except for the very critical regime, also for 
the sc(110) surface (Fig. 4). Furthermore, we define the 
ratio r 2 — z$/z. Two-dimensional Fourier transforma- 
tion of the renormalized hopping matrix T = \/ZT\/Z 
yields a tridiagonal matrix T aa > (k) in the layer indices 
a, a' — 1, ...,N±. For a, a' > 2 the non-zero elements 
are given by T aa (k) = ze||(k) and T act±1 (k) = ze±(k). 
For the sc(100) surface the parallel and perpendicular 
dispersions read pq ]: ey(k) = 2t(cos(k x a) + cos(k y a)) 
and e_i_(k) = t. At the very surface we have: Tn(k) = 
r 2 ze||(k) and T 12 / 2 i(k) = rze±(k). The specific sur- 
face renormalization (r 2 < 1) leads to two (localization) 
effects relevant on the energy scale zW: a diminished 
parallel dispersion in the top layer and a tendency to 
decouple the top-layer from the bulk spectrum. 

As a consequence, surface excitations split off the con- 
tinuum of bulk excitations in the SBZ: For tridiagonal 
T aa i (k) the k-resolved surface Green function is readily 
calculated 1241; we have: 



G (k,£) 



- a T VP ~ a Y ~ 4 & 2 ) 

(for ± Re(E - a) > 0) , (8) 



with a = zeii(k), b = ze^(k) for the unperturbed surface 
(r 2 = 1). Thus, for r 2 = 1 the coherent part of the 
surface spectral density — ImGo(k, E + iQ + )/ir is semi- 
elliptical for each k with band edges at E = a ± 2b. 
The surface Green function G for the perturbed surface 



< 1) can be expressed in terms of Go as 



G(k,E) 



zs 



E 



2 & 2 G (k,£) 



(9) 



A surface excitation is characterized by an energy that 
(for a given wave vector k of the SBZ) falls outside the 
continuum of bulk excitations described by Go- A con- 
dition for the existence of surface excitations is there- 
fore given by \E$ — a\ > 2b where Eg is the excitation 
energy obtained from the perturbed surface Green func- 
tion: G(k,£ = ^s)- 1 = 0. With G (k,E) = ±z/b at 
the band edges E = a ± 2b this is equivalent to: 



> 2b 2 ± b(a - <z ) - b 2 



(10) 



where ao 



-II ( k ) 



2 a and bo = rzej_(k) = rb. Eq. 



(10) generally tells us whether or not a discrete eigen- 
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Fig. 5. Upper left: z(U), zs(U) and the ratio r 2 
for the sc(100) surface. Upper right: (2 — r 2 )/(l — 
function of r 2 . Lower right: Ratio between parallel and per- 
pendicular dispersion along high-symmetry directions in the 
two-dimensional surface Brillouin zone (SBZ). At k points 
with (2 — r 2 )/(l — r 2 ) < jey (k)/ex(k)| surface excitations 
are split off the bulk continuum. Lower left: Coherent part 
of the quasi-particle spectrum projected onto the SBZ cal- 
culated from Eqs. (^) and (|^). Filled area: bulk continuum. 
Solid line: surface excitation for r 2 = 0.25 corresponding to 
U = Up- 



value is split off the eigenvalue continuum of a semi- 
infinite tridiagonal matrix with modified "surface" pa- 
rameters. For the present case it immediately leads to 
the following simple criterion for the existence of a sur- 
face excitation: 



1 



< 



|(k) 



(11) 



Hereby, the ratio r 2 = z$/z is related to the ratio be- 
tween the free dispersions. 

Fig. 5 visualizes the criterion for the case of the 
sc(100) surface. At small U , the renormalization factor 
r 2 is too close to 1 to generate surface excitations. These 
appear above U ss lit and split off the bulk continuum 



at the r and M high-symmetry points in the SBZ. The 
minimum interaction strength is well below U C 2- Finite- 
size effects are of no importance here. For U i— > U C 2 
surface excitations are found in a fairly extended region 
around T and M . In the insulating phase the whole co- 
herent spectrum diappears. The sc(110) surface is even 
more favorable for surface excitations. Eq. (|ll|) always 
holds along the k y = n/^/2a direction where ej_(k) = 
pr| . For the sc(lll) surface z a deviates from z also in 
some subsurface layers. Calculating the coherent surface 
spectrum from the DMFT results for z a (U) shows that 
surface excitations do not exist for any U . 

As can be seen in Fig. 5, the surface excitation splits 
off the bulk continuum away from the Fermi energy. Due 
to quasi-particle damping the surface peak in the spec- 
trum will therefore aquire a finite width which is of the 
order AE w z$ ■ |7s|-Eg where 73 is the damping coef- 
ficient at the surface defined in Eq. (|J). Since we have 
E§ w zgi for the position of the excitation, the width 
is AE w z|| 7s |t 2 . For U i-> U c2 {z ^ 0) the E 2 co- 
efficient in the expansion of the self-energy diverges as 
7 ~ — l/z 2 W as has been shown in Ref. (r| for the 
d = 00 Bethe lattice. This result may be used to esti- 
mate AE < zgt. Comparing AE with the energetic dis- 
tance between the surface peak and the edge of the bulk 
continuum, |J5g — £b| ~ zt (Fig. 5), shows both quanti- 
ties to be of the same order of magnitude for U i—> U C 2- 
We conclude that the surface peak cannot be fully sep- 
arated in energy from the bulk features in the surface- 
projected spectrum. 

A complete energetic separation between bulk and 
surface excitations is possible if E$ crosses the Fermi 
energy as a function of k. The dispersion of the surface 
mode can be calculated from (H) and (0) to be: 



E s (k) = i 2 r 2 ej|(k)±zr 2 i 



h(k) 2 



-(k) 2 



1 



(12) 



Setting Es(k) = yields ej_(k) = 0. As mentioned 
above, this is fulfilled e. g. along the k y = n/\/2a 
direction for the sc(110) surface. 



VI. CONCLUSION 

While the investigation of the surface electronic struc- 
ture of a metal has a long tradition, there are only 
few attempts to account for correlation effects beyond 
the Hartree-Fock approximation or the (ab initio) local- 
density approach (see |^0[ and references therein). 
The present paper has focussed on the Hubbard model 
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as a standard model for correlated itinerant electrons 
and on the Mott transition as a genuine correlation ef- 
fect. Adapting the dynamical mean-field theory for gen- 
eral film geometries and using the exact-diagonalization 
approach, we have investigated the specific surface prop- 
erties of the T — semi-infinite Hubbard model at half- 
filling. 

At the surface electron-correlation effects are found 
to be significantly more pronounced compared with the 
bulk. The top-layer quasi-particle weight is consider- 
ably reduced for all surfaces and interaction strengths. 
Surface correlation effects are the stronger the larger is 
the reduction of the surface coordination number. For 
the comparatively open sc(lll) surface there is a pro- 
nounced layer dependence of the quasi-particle weight 
z a for weak and intermediate couplings. As U ap- 
proaches the critical region, the layer dependence be- 
comes monotonous in all cases. 

The reduced surface coordination number tends to 
drive the system to a phase with an insulating sur- 
face on top of a metallic bulk. Eventually, however, 
the low-energy bulk excitations induce a finite spectral 
weight at the Fermi energy in the top-layer density of 
states. Consequently, all layer-dependent quasi-particle 
weights vanish at the same coupling strength U C 2] there 
is a unique transition point. This excludes a modified 
surface critical interaction and thus the existence of a 
surface phase. 

On the other hand, a different localization effect has 
been found: There are special regions in the surface Bril- 
louin zone, where the low-energy excitations in the top 
layer cannot couple to the bulk modes and are confined 
to two-dimensional lateral propagation. These surface 
excitations are of particular interest since their nature is 
rather different compared with the well-known Tamm- 
and Shockley-type surface states (£8), ^9|. Opposed to 
the latter they are caused by electron correlations ex- 
clusively. Quite generally, as long as the concept of a 
highly renormalized Fermi liquid applies, the specific 
surface renormalization of the effective mass tends to 
generate excitations split off the bulk continuum. This 
mechanism is sufficiently general to survive also at finite 
temperatures and non-bipartite lattices where antifcr- 
romagnetic order is expected to be suppressed. High- 
resolution photoemission from single-crystal samples of 
transition-metal oxides should be appropriate to detect 
the excitation by suitably tuning the Mott transition. 

The present study has been restricted to surface 
geometries and uniform model parameters. It is an 
open question whether or not there is a surface phase 
if the hopping or the interaction strength at the very 



surface is modified. Possibly, a metallic surface coex- 
isting with a Mott-insulating bulk can be found for a 
strongly reduced surface U. Future investigations may 
also concern thin-film geometries where the study of 
the thickness dependence of the critical interaction is of 
particular interest. 
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